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Abstract 

In this letter, we report the existence of a novel type of explode-decay dromions, 
which are exponentially localized coherent structures whose amplitude varies with 
time, through Hirota method for a nonisospectral Davey-Stewartson equation I dis- 
cussed recently by Jiang. Using suitable transformations, we also point out such 
solutions also exist for the isospectral Davey-Stewartson I equation itself for a careful 
choice of the potentials. 

The identification of exponentially localized structures namely the so called "dromions" 
in the Davey-Stewartson I (DSI) equation [1] has given a new direction to the study of non- 
linear partial differential equations (pdes) in (2+1) dimensions. The subsequent investiga- 
tions of other physically and mathematically important (2+1) dimensional nonlinear pdes 
[2-10] have enriched the family of integrable nonlinear pdes, exhibiting considerable rich- 
ness in the structure of the solutions. However, very little is known about non-isospectral 
(2+1) dimensional nonlinear pdes in this regard. Now one may ask, do exponentially lo- 
calized structures exist in the later systems also, and if so, what are their characteristics? 
This letter is a modest attempt in this direction in which the existence of such localized 
solutions is pointed out for DSI equation with nonuniform or inhomogeneous terms. 

In this letter, we consider a non-isospectral DSI equation introduced by Jiang [11] of 
the form 



ilt + ^(Qxx + Qyy) + 



2 ^ wjv 2 
i{uiy + a\)q y - iwi(xq) x + 2(uj x + a )q = 0, 



(la) 



• -If 

t^t „ y"xx ' fyy 



+u 2 {ri,t) - ^<9 ? 1 (qr) v -d r) 1 {qr)t: 



(lb) 

i(u>iy + ai)r y - iuJi(xr) x + 2(uj x + a )r = 0, £ = x + y, n = x-y, 
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where too, wi, ao and a\ are all real constants and u±(^,t) and 112(1], t) are the so called 
boundaries (arbitrary functions). The above set of equations has been investigated via 
inverse scattering transform method [11] and an integro-differential equation for the time 
evolution of the scattering data by virtue of the time dependence of the scattering param- 
eters has been brought out. Under the reduction r = — q*, equation (1) reduces to the 
following form 



kt + Qtf + <lrm + 



(2) 



-zwi[£<?£ + r]q v ] - iai[q^ - q v ] + [2a + wo(£ + v) ~ iuJ i]q = 0. 



In equation (2), the parameter ao can be scaled away by introducing the transformation 
q = qexp (—2iaot). Then, the above equation can be equivalently written as 

\ (I 1 \ 2 ) r , = (3a) 

\ (I Q 1% = V v , (3b) 
iqt + q^ + q m + (U + V)q-iuJi(^ + r]q v )-ia 1 (q^-q 71 ) + (w (£ + v) ~ iuJ i)q = 0.(4) 

As the complete integrability parallelling to that of dynamical systems under isospectral 
flows is not obvious for the nonisospectral problems in (2+1) dimensions, we address our- 
selves only to the nature of the solutions of the above equation rather than its integrability 
property. For this purpose, we bilinearize equation (3) and obtain the solutions using the 
Hirota method. 

To bilinearise the above equation, we effect the following dependent variable transfor- 
mation 

q = j, U = 2d m logF, F = 2%logF, (5) 
so that equations (3) and (4) get converted into the following Hirota form, 
iDt + DZ + Dl-iuiiSDz + riDr,) 

1 (6a) 
-iai(D € - £>„) + (wo[£ + 77] - iui)] G ■ F = 0, 

2D^D V F-F=\ G | 2 . (6b) 

We now introduce the following power series expansion 

G = eg^ + e 3 5 ( 3 ) + . . . , (7a) 

F = 1 + e 2 /^ + e 4 / {4) + • • • , (7b) 

into the bilinear form (6), where e is a small parameter. Collecting the various powers 
of e, we get the following set of equations, 

0(e): i(5 (1) ) t + (5 (1) )^ + (5 (1) ) w -^ 1 [^ (1) ) $ + r?(5 (1) ),] 

(8a) 

-i ai [(gW)t - (g(%] + (uo(t + r,) - iuJgW = 0, 
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0(e 2 ): Mf^ = g^9 {1) *, (8b) 



0(e 3 



iD t + D| + £)2 - iwx^Df + r,D v ) 

(8c) 



-iai{D^ - D v ) + (w [£ + ??] - iui) 



+ ( 5 (3) +5 (l). / (2) )=0) 



0( e 4 ) : 2£> € A,(2/< 4 > + /( 2 >./ (2) ) = ^ty 1 )* + ^ V 3) \ (8d) 



and so on. Solving (8a), we obtain the simplest "plane wave" solution 

AT 

5 (i) = e* , = % + «j (t)»7 + / fij 

7=1 7 



(9) 



where the spectral parameters kj(t) and Zj(t) evolve in an identical fashion obeying the 
time evolution equation 

i(kj) t - i(ui)kj = i(lj) t - = -vo, (10) 

and the "dispersion" relation is given by 

iQj(t) + (kjit)) 2 + (lj(t)) 2 - ia 1 {k j (t) - lj(t)) - iui = 0. (11) 
To construct one soliton solution, we take N = 1 so that we have 

gW= e x\ (12) 
and hence the solution of (8b) becomes 

; ' " l6k 1R (t)hR(ty 1 j 

Substituting and /( 2 ) in equations (8c) and (8d), one can show that = for j > 3 
and = for j > 4 provided the spectral parameters obey the following equation 

{k\ R h R )t = 2uikiRhR. (14) 

Considering the above equation alongwith equation (10), separating out the real and 
imaginary parts of the later, the time evolution of the spectral parameters can be obtained 

as 

fcifl(t) = k 1R {Q)e^\ l 1R (t) = *i fl (0)e wi * (15a) 

and 

ku(t) = kuiOje^ - ^, l u (t) = /i/(0)e^* - ^, (15b) 

where ki R (0), ku(0), Iir(0) and lu(0) are all constants. Using (5), (7), (12) and (13), the 
physical field variable q of equation (4) is driven by the envelope soliton (line soliton) 



q = 2^/k lR {t)l lR {t) sech ( X iR + V>K Xli '. (16) 
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It is interesting to note that the amplitude of the above soliton solution varies with time 
by virtue of equation (15). Similarly, the potentials U and V are driven by the line soliton 
solutions whose amplitude again varies with time as 

U = (k 1R (t)f sech 2 ( X ifl + V), (17a) 
V = (l 1R (t)) 2 sech 2 (xiR + 4>)- (17b) 

As it is clear from equation (3) that the boundaries are specified by the arbitrary 
functions 112(1], t) and u\(£,t) which drive the potentials U and V even in the absence 
of the field q, as in the case of isospectral (uniform) DSI equation, one can expect the 
presence of exponentially localized solutions in the system provided one properly takes care 
of the time evolution of the spectral parameters. To generate a (1,1) dromion solution, we 
take the ansatz 



(18) 



X2 = hv + J £l 2 (t)dt, Qi(t) = ikf + aiki, Q 2 {t) = ill — aih. 

where a, j3, 7 and 5 are parameters. Substituting (18) into (6b), we obtain 

G = pe Xl+X2 , p is real, (19a) 

p 2 = 8k 1R (t)l 1R (t)(5j -af3). (19b) 



Hence, the dromion solution now becomes 

pe Xi+X2 



§ _l_ ae Xl+x i + (3e X2+x 2 + ^ e xi+xl+X2+xl ' 



(20) 



It can be easily observed from the above solution that the amplitude of the dromion 
solution has to evolve in time obeying the time evolution equation 

p t - uJip = 0, (21) 

so that solution (20) satisfies equation (6). The above equation has the solution 

p = poe^, (22) 

where po is a constant. Thus, the amplitude of the dromion solution varies with time in 
the nonisospectral DSI equation unlike the isospectral DSI equation where it remains a 
constant. It should be mentioned that both the spectral parameters kj and lj as well as 
the amplitude of the dromion solution are governed by the same time evolution equation 
(differing only in integration constants). Thus, depending on the nature of the parameters 
involved, the amplitude of the dromion solution (20) either grows (explodes) or decays 
with time. A typical example is shown in Figs. 1-3. We call these types of solutions 
as "explode-decay dromions" which is reminiscent of the explode-decay solitons of the 
inhomogeneous (1+1) dimensional nonlinear pdes [12]. To our knowledge this seems to 
be the first instance of such a localized solution being realized in a (2+1) dimensional 
nonlinear pde. 
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Finally, we wish to point out that the above type of localized solutions with time 
varying amplitude do exist for the isospectral Davey-Stewartson I equation also. Using 
the following transformations 



q = q exp 



U 



1 



-aiui - u 77, 



V = V + -uli 2 - 
4 1 



the nonisospectral DSI equation (3)-(4) reduces to the isospectral DSI equation 
m + q# + q m + (u + V) q = 0, 



1 I - |2 



q k 



(23a) 
(23b) 
(23c) 

(24a) 
(24b) 
(24c) 



It is well known that equation (24) admits exponentially localized solutions with constant 
amplitude by driving U and V by sech 2 potentials [1, 8]. However, by virtue of the above 
transformation (23), it is now evident that such localized solutions with time varying 
amplitude do exist for the isospectral DSI equation also for a careful choice of the potentials 
indicated by the transformations (23b) and (23c). However, this choice of the potential is 
not very obvious but for the above nonisospectral case by virtue of the transformations (23) 
and hence such localized solutions whose amplitude either grows or decays with time have 
eluded earlier observation. 

In this letter, we have generated a new class of localized coherent structures to DSI 
equation with inhomogeneous terms known as "explode-decay dromions" whose amplitude 
varies with time unlike the basic dromions. It remains to be seen how the multiexplode- 
decay dromions would interact in the context of the variation of the spectral parameter and 
this remains as an open question. We have also indicated the possibility of the existence 
of such solutions for the isospectral DSI equation itself. 
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Figure Captions 

Fig. 1: The time evolution of explode-decay dromion at t = —0.5 
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Fig. 2: | q | at t = 



Fig. 3: | q \ at t = 0.2 
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